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Abstract 



An approximation method, namely, the Extended Wronskian Determinant Ap- 
proach, is suggested to study the one-dimensional Dirac equation. An integral equation 
which can be solved by iterative procedure to find the wave functions is established. 
We employ this approach to study the one-dimensional Dirac equation with one-well 
potential, and give the energy levels and wave functions up to the first order itera- 
tive approximation. For double-well potential, the energy levels up to the first order 
approximation are given. 
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1 Introduction 



Since the majority of problems in quantum mechanics lead to very complex equa- 
tions and cannot be solved exactly, the approximation methods for solving Schrodinger 
equation are very important for practical use. Recently, Friedberg, Lee and Zhao(FLZ) 
suggested a new method to decide the wave functions and energy levels of Schrodinger 
equation by quadrature along a single trajectory [1-3]. Giving a good trial wave function 
^ for one-dimensional(lD) Schrodinger equation with potential V, 

1 d^^ 

they considered a perturbation equation 

Id^^ 



2 dx2 



+ (y + w)^ = {E + e)^ (2) 



where w is the perturbation potential and e the energy shift. By means of the Green's 
function, they constructed an integral equation and solved it by iterative procedure to 
find the true ground wave function and the corresponding energy level. The integral 
equation for wave function ^ reads 

f+oo r+OD 

$ = $-2$/ ^''^{y)dy / <S>{z){w - z)'^{z)dz (3) 

J X J y 

and the corresponding energy is 

■-co 

— (4) 







^^dx 







Employing a double-well potential^ and other potentials[4, 5] they gave the ground 
state energy up to the first order approximation of iterative procedure. 

On the other hand, as is well known, the behavior of ID Schrodinger equation is 
determined completely by its corresponding Wronskian determinant [6]. Therefore, a 
Wronskian Determinant Approach(WDA) is suggested by us [7] to study the energy 
levels and the wave functions of ID Schrodinger equation. We established an integral 
equation and employed this equation to study the same problem of double-well poten- 
tial. A series of expansion of ground state energy up to the second order approximation 
of iterative procedure is given. We found that the series of energy given by FLZ method 
and the WDA up to the order of 0{g^^) where g is the strong coupling constant are 
the same. This result confirms that the WDA is a successful method to investigate the 
ID Schrodinger equation. 

The objective of this paper is to extend the WDA to study the ID Dirac equation. 
Based on general considerations and an extended Wronskian determinant, we establish 
an integral equation for Dirac wave function and give the "resonant condition" to deter- 
mine the energy level. Our approach may be considered as a relativistic generalization 
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of the WDA. We will show the formulae of our Extended Wronskian Determinant Ap- 

proach(EWDA) in next section. In section III, as an example, we will use EWDA to 
study the ID Dirac equation with perturbative one-well potential. We will solve the 
integral equation by iterative procedure up to the first order to find the energy shifts 
and the wave functions. Since the double-well potential for Schrodinger equation and, 
in particular, for Dirac equation has wide applications in physics, for example, soliton 
bag model, instantons and anti-instantons and etc[8-12], we will address the energy 
levels and the wave functions for double-well potential in section IV. The last section 
is a summary. 
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2 The Extended Wronskian Determinant Approach 



The ID Dirac equation reads 

{ap + Pm)ilj{y) = {E-V)ilj{y) (5) 

where if) = (^) is a 2-component spinor. Choose a = (T2, $ = ^3 and define 

. _ / m + E-V\ 

^=[m-E + V J 

then Eq. © could be rewritten as 

ij' (y) - AV (y) = (7) 

Eq. ((T)) is called the homogeneous equation in our approach, We need its specific solution 
as a zeroth order approximation to construct the solution of the perturbation equation. 
Assume a perturbative potential w is added to Eq. © 

V = V + w (8) 

and the eigenenergy is shifted to 

£ = E + e (9) 

We get the perturbation equation 

V^' (y) -Ai^ = {e-w) U (10) 

where 

;) (") 

Eq. (|1U|) is called an inhomogeneous equation and the inhomogeneous term is {e — w) iip 
in our approach. 

As is well known, in ordinary differential equation theory, with a specific solu- 
tion of a homogeneous differential equation given, we can try to find another linearly 
independent solution of this equation and the general solutions of the corresponding 
inhomogeneous differential equation through Wronskian determinant calculation. This 
procedure can also be employed to discuss the solution of ID Dirac equation. 

Suppose ■01 = ("^) is a specific solution of Eq. Q. It is a normalized bound state 
and satisfies the boundary condition 

lim ipi = 0. 

y— »±oo 

Now we construct another solution ip2 = ("j) which is a linearly independent solution 
with tpi by using the extended Wronskian determinant approach. Obviously we have 

{y2 -U2) ip'i - (f 1 -ui) ip2 = {v2 -U2) Atpi - [vi -ui) A'ip2 (12) 
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A straightforward calculation can show the right hand side of Eq. H12|l be zero, we get 

{uiV2 - U2V1)' = (13) 



Define a matrix U as 



We find 



U 



Ul Ui 
Vl V2 



detU = C (Const.) 



(14) 



(15) 



Here detU is called an extended Wronskian determinant. The constant C is nonzero 
because ipi and ■02 are linearly independent. Since "01 satisfies the boundary condition 

lim ipi = 0, 

y— »±oo 

we get that ■02 niust be divergent when y —i- ±00 because of Eq. ((T3|). Eqs. (|T^ and 
yield 

u,V2 - ^_l^y v'2n = C (16) 

Assume V2{y) has the form V2{y) = h{y)vi{y) where h{y) is an unknown function 
which can be solved from Eqs. (|16|) and 0, we find 

m- E + V 



h{y) = -C dy 



and 



V2 



-Cvi / dy 



m- E + V 



(17) 
(18) 



Similar procedure can be used to construct 1*2, we finally obtain 



1p2 



-c 



Ul I dy- 



Vl / dy 



-m- E + V 
+m - £' + y 



(19) 
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■02 dose not satisfy the boundary condition of V'l •, because of the uniqueness theorem. 

Now we are in a position to find the solution = (^) of the corresponding inho- 
mogeneous perturbation equation ((TU)) . Suppose ^0 has the form 



= (J)) = U(y)<7(y) 



(20) 



where g[y) is a column matrix. Substituting Eqs. (|2Uj) and (|14j) into Eq. (|1U|) we obtain 

g' = {e- u;)U-ifU5 = (21) 
where is the inverse matrix of U. A simple calculation can show that 

'UIU2 + V1V2 



1 



B (e - u;, y) = (e - ?i;)U"^fU = — (e-w 



-(nf + 't;f) -(uiU2^ V1V2)/ 



(22) 
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and the wave function 

i) = t j AyBg = t j dy{e - w)ij-^i^ (23) 

or, explicitly, 

Eq. H24() is an integral equation and can be solved by iterative procedure step by step. It 
can be proven that the choice of the constant C does not affect on our result remarkably. 
For simplicity, we choose C = 1 in the following calculations. The integration bound of 
Eq. (|24j) is determined by the boundary condition and the normalized condition. Based 
on Eq. H21|) . we obtain an iterative sequence of g as 

{e-w,y) = g{e-w,a) 

5f(^) {e - w,y) = g{e - w,a) + / dy'B (e - w, y') g^^'^ (e - w, y') 

J a 



(25) 



ry 

5^") {e-w,y)=g{e-w,a)+ / dy'B (e - w, y') g^"^'^^ {e - w, y') 

J a 

The corresponding formal solution of wave function tp is given by Eq. H20|l . 
Now we turn to discuss the energy shift e. We get 

^\ = ( ^^^^ ~^ ^2^2 
v) \vigi+V2g2 



(26) 



from Eq. (|2()j) , where g = ( ) . Noting that the boundary conditions for and ijj2 

\92; 



are 



and 



lim ipi = lim ( 



lim il:2 = lim = oo 

y— >itoo y— +±00 \ V2 



respectively, we have 



lirn 52 = (27) 

y— >±oo 



According to Eqs. (jH]) and ((221) we find 

1 fy 



g2{e -w,y) = g2{e - w,a) + j dy {e - w)[-{u\ + vl)gi - {uiU2 + viV2)g2] (28) 
By means of Eq. p7|l , we get 

dy (e - w) [{u\ + vf) gi + {uiU2 + ^1^2) 52] = (29) 



/ 

J — ( 
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Eq. H29|l is called the "resonant condition". The energy shift reads 



+00 



dyw 



{uj + vf) gi + {uiU2 + V1V2) 92 



dy 



{uj + vf) gi + {uiU2 + V1V2) 92 



(30) 



Eqs. (|24jl . (|25jl . (|3U|) as well as the normalization condition of V' are the basic formulae 
in EWDA. The zeroth order wave function can be taken as V'^''^ = V'l when w is a 
perturbation potential. Starting from we can find the wave function and the 

energy level up to any order which we need in principle by iterative procedure. 
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3 One- Well Potential 



3.1 Square- Well Potential 

As an application and example, now we study how to find the bound state wave func- 
tions and the energy levels of ID Dirac equation with one well-potential. As a homo- 
geneous equation, we study the square-well potential first. 



The ID Dirac equation with square-well potential 

V{y) = -Vo9{y + ^)e[-y + ^) 
where 6 (y) is the step function 

\0 (!,<0) 

(Fig. 1) can be solved exactly [13]. The bound state solution is 



Ui 



K / 



+ 



+ 



COS py 

m-E-Vn 



sm py J 



1 

m-E 

K / 



UB + 



sm py 
m-E-Vo 



uc 



cos py. 



in which 



K = a/ m? — E'^ 



(— m < E < m) 



p = \J{E + V^f - m? (m - Vb < E) 
The coefHcients ua,ub,uc and ud m one case are given by 



UA = ud = UBe^ cos -pa 



UB 



m 1 



m + E K 



(1 + cos pa) 



E + Vo 



-a + 



m 



1 



E + Vo + m E + Vo + mp 
uc = 



sm pa 



and the energy spectrum satisfies 



cot -pa 
2 



(m + E){E + Vo- m) 
{m -E){E + Vo + m) 



(31) 



(32) 



O(y + ^]0i-y + 



(33) 
(34) 

(35) 

(36) 

(37) 
(38) 

(39) 
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The upper component ui and the lower component vi of this solution -01^/ are 
even symmetric and odd symmetric respectively. In another case, the coefficients of 
the alternative solution tpi jj in which the symmetry exchanges are given by 



ua = —U£) = —uce^ sm -pa 



(40) 



uc 



m 



UB = 



^1 — cospaj + — — a 



m 



1 



m + E K E + Vq + m E + VQ + mp 

and the energy spectrum satisfies 



sm. pa 



(41) 
(42) 



tan -pa 



{m + E){E + Vo- m) 
(m -E)(E + Vo + m) 



(43) 



Another linearly independent solution 1^2 can be obtained by using Eq. H19|) . The 
result reads 



U2 



+ 



^e~^ye i-y - ^ 



' m+E+Vo 



sm py\ 1 



m+E+Vo 



COS py 



UB 



+ 



p 

sin py 



cos py\ 1 
/ uc 



-^+2)^^+2 



2^ I J UD 2 



(44) 

where ua,ub,uc and ud are given by Eqs. (38-40) for ^jJlJ and Eqs. (42-44) for 'ipiji 
respectively. 



3.2 Perturbation of the Square Well Potential 

Now we go to discuss the inhomogeneous equation (|10|) . The one- well potential is given 

by 

>{y^l)o{-v^l) (45) 

(Fig. 2) and the perturbation potential reads 

-(.) = ^y^^(y + |)^(-y+^) (46) 

We calculate the first order energy shift £^^\ Using ipi as the zeroth order non- 
perturbation wave function V'^"^ then 5(0) = Q from Eq. Substituting Eq. ^ 



V = -Vq 
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and g^'^^ into Eq. H3U|) . we obtain 



4Vb 

2 



-(1) _ « 



m ( a a 1 

„ - , „ — sm va H ^ cos pa 5- sm pa 

£; + Vb + m 12 E + VQ + m\Ap ^ 2p^ ^ 2p^ ^ 



E + m K 



for ip^^^ = ipi J, and 



ml, , / ^ + Vb 

-— (1 + cos pa) + ( ^ , — : a + 



m 



1 



„ - ' ~r7 smpa 



(47) 



4Vb 



-(1) 



£^ + ^0 1 



m 



^ + Vo + m 12 ^ + Vb + m 



a 



I — sin pa + —J cos pa — — ^ smpa 



1 
2p 



m 1 ^ , ^ + ^0 

(1 — cos pa) + I — — a 



E + m K 



m 



E + VQ + m E + Vo + mp 



sin pa 



(48) 



for = il^i^ II. 

The wave function for one-well potential has been calculated up to the first order 
iterative approximation. The results are 



(1) 



(1) 



/ 1 / ^ m 
I -i-Ey + — 
K \ 2k 

1 / m 
~m + E\ ^^2k. 



-y - 



+ 



1 
p 



e*-^-* {A cospy + C sinpy) ^ {B cospy + D sinpy) 



4Fo 



\E + Vo + 



m 



e^^^ {—Asinpy + C cospy) ^ {—B sinpy + cospy) 



/ 1 / ^ m 



for V''''''^ = V'l, /) a.nd 




1 / p 



(49) 



-y - 



+ 



e*-^-* (A cospy— C sin py) ^ {B cospy— D sin py) 



1 



e^^^ (+^sin py + C cospy) ^ (+i? sin py + cospy) 



uc- 



y + 



/ 1 

V 



^ m 
Ey + — 
K \ 2k 



1 / 7-1 

m + E\ ^^2k, 



(50) 
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for tp^^^ = -01^ II where 



m 



A = —A— cos 2py 
2p 



B = -X( —y cos 2py - sin 2py 



D = 



2p' 
C=- 



cos 2py 



m 



{E + Vo)y + \—sm2py 
2p 



m 



^ {E + Vo)y'^ + X [ ^y^ sin 2py + -^y cos 2py - sin 2py 



m 



m 



and 



A = 



f + 1 (for V'^") = Vi, /) 



4pd 



[ - 1 (for = V'l, //) 
C^^) is determined by normalization condition. 



(51) 
(52) 
(53) 
(54) 

(55) 
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4 Double- Well Potential 



The ID Schrodingcr equation with double-well potential has been investigated by FLZ 
method [3] , WDA and variational method [7] previously. We now apply the EWDA to 
the double-well potential of the ID Dirac equation and give the energy eigenvalues to 
the first order approximation. 



4.1 Double- Well and Double-Square- Well Potentials 



Given a symmetrical double-well potential 

2 



(Fig. 3) with the minimum of V being — /x, define 



then we have 

With r and rj, instead of r) and ^, V can be expressed as 



(56) 



(57) 
(58) 

(59) 



An exactly solvable doublc-square-well potential, whose bound states can serve as 
the zeroth order wave function for iteration procedure, is chosen in the homogeneous 
equation as. 



in which 



V =Vo9 i-y -a)-^9{y + a)9 {-y - b) + i^O {y + b) 6 {-y + b) 
- (10 {y - b)e i-y + a) + Voe (y - a) 



Vb = V(-a) = V(a) 



(60) 



(61) 



where a and 6(0 < b < a) are parameters which can be adjusted to minimize the energy 
shift e^^l 



The wave function of a bound state solution to the Dirac equation with the double- 
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square-well potential is 



Ui 



where 



+ I m-E + Vo j UAe^^e {-y - a) 

ub cos Ay + ub sin Xy 
+ I m — E — fi 



A 



{ub sin Xy — Ub cos Xy) 



+ 



cosh qy 
m — E + u 



'iy + b)9 i-y - b) 



sinh qy 



uc+\m-E + v ^ \uc 
smh qy j \ cosh qy 



{y + h)9 {-y 



Ub cos Xy + ub sin Xy 
+ I m — E — fj, 



X 



{ub sin Xy — Ub cos Ay) 



'{y-b)e{-y + a) 



1 



+ rn-E + Vo]'^Ee''ye{y- a) 



-{E- Vo)^ {-m -Vo<E <m + Vo) 



(62) 
(63) 



A = ■\/ —m^ + (£■ + /x)2 (£^ < — m — jj, or m — n < E) 



(64) 



g = \/ — [E — z/)2 (— m — v<E<m + v) 
If Eq. ()65p dose not satisfy when (i? — v)"^ > m-^, we define 



p = \J —m^ + (S — < — IJ, or m + < E) 



(65) 



(66) 



By replacing g by ip and -uc" by —iuc the wave function ()62() and all the following 
expressions maintain their forms. 

Similar to the square-well case, the coefficients ua,ub,ub,uc,uc,ub,ub and ue 
of the solution j with even symmetric upper component ui and odd symmetric 
lower component vi are given by 



Ub = Ub = Oiuc 



(67) 



Ub = —Ub = —auc 



UA = UE = juc 



(68) 
(69) 
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Uc 



Am 



m + E-VQ 
4m 



4m 



m + E + /Lt^ ^ m + £; + Ai^ ^2^ ' 



2m 1 , , 2{E-v) ' 

-aaL H —smhzqb H b 

m + E — u 



m + E + 11 m + E — u2q 

and the energy spectrum satisfies 

K m — E — jjL 
tanhg6^ ^-E + V, 



X 



tan X{b — a) — 1 



A 



m — E — fx 



tan A (6 — o) + 



m - + Vb 



m — E + n 



(70) 



(71) 



Correspondingly, the coefficients of tlie alternative solution ij^i^ u with the ex- 
changed symmetry of components are given by 



UB = -UD = -fiuc 



(72) 



UB =UD = fiuc 



UA = -UE = -Sue 



(73) 
(74) 



uc 



4m ^9 ^ 4m 



m + E -Vq m + E + jjL 



m + E + n 



nnr ^m 1 . , „ , 2 {E - v) ' 

PPL + sinh 2qb ^— — -b 

... ,771 _.o„ 



m + E + ji m + E — v2q 

and the energy spectrum satisfies 

K m — E — jj, 

coth qb = 



m- E + Vq a 



tan A (6 — a) — 1 



A 



tan A (6 — a) + 



m — E + fi 



— E — ji m — E + Vo^ 

In the above expressions we have defined 

a = cosh qb cos Xb -\ ^ — — " sinh qb sin Xb 



a = cosh qb sin Xb — 
P = sinh qb cos A6 + 



A 


m 


-E + v 


m — E 




q 


X 


m 


-E + v 


m — E 




q 


X 


rn 


-~E + u 



m — E — ji 



sinh qb cos Xb 
cosh Q't sin A6 



(75) 



(76) 



(77) 

(78) 
(79) 
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7 = e 



5 = 6" 



(3 = sinh qh sin Xh 



cosh qb cos \{b — a) + 



sinh g6 cos A(6 — a) + 



m — E + V 



m 


-E 


-/^ 


Q 




A 


m 


- E + v 


m 




- M 


q 




A 


m 


-E + v 


m 


- E 




q 



cosh g6 cos A6 (80) 

sinh g6 cos A (6 — a) (81) 
cosh qb sin A(6 — a) (82) 



and C, and L given later in Eqs. (90), (93) and (96). 



Another hnearly independent solution, though it is not necessary for deriving the 
first order energy shift. But it is involved in higher order calculations. The result, 
which can be constructed by Eq. (fT!l|) . reads 







C) 









-m- E + Vo'^ 

K 
1 



1 

( 

UA 



i-y-a) 



' m + E + fi 



+ 



+ 



+ 



A 



{ub sin Xy — ub cos Ay) 



{ub cos Xy + Ub sin Xy) 



9{y + a)0i-y-b) 



'm + E-v . 1 \ / m + E-v , , ^ 
smh I i / cosh qy \ i 

q — ^ q \ — 

cosh qy / \ ~ sinh qy ' 



^{y + b)6{-y + b) 



' m + E + n 



X 



{ud sm Xy — U£) cos Xy)\ 1 



+ 



UD cos Xy + UD sin Xy 

e'^ye {y - a) 



iy-b)e {-y + a) 



^ ' m + E-Vp . ^ 



K 
1 



UE 



(83) 

where u^, ub,ub,uc, uq, ud,ud and are given by Eqs. (69-72) for ip2, i and Eqs. (74- 
77) for ■ip2^ II respectively. 



4.2 First Order Energy Shift for the Double- Well Potential 

The perturbation potential on the double-square-well potential reads 

w{y) = Viv) - V{y) (84) 

Using tpi as the zeroth order tp^^^ in iterative procedure, and g^^^ = (q) from 
Eq. (|2U|) . the first order energy shift e^^^ is obtained by substituting Eq. ()62() and g^'^^ 
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into Eq. (jHOJ, we get 

4m 



+ 



+ 



+ 



+ 



m + E -Vq 

Am 
m + E + fj, 

ijE + fi) 
m + E + fj, 

Am 



A 



D 



V 1 



1 \ 1/ 1 



(a^ _ _ 1 + (a^ - 6^) + + ^) (a - 6) 



^2 y ^2 g 



m + E + n 



UBUB 



^2^4 



J 



1 + ^ ) 



+ 



+ 



2m 



m + — 
2{E-v) ,^,2 



M 



m + E — V 
where 



uc - uc) 



V 1, 



^2^45 



1 + ^)4^ 

r^y r/^ 



1 ^ ^ ^ ^5^' 

j-2 j ^23 



.a a 3a 3a 3 \ .2^0 



(85) 
(86) 



^ 1 /a^ a 1 
S = :t — + ^ + 



2 V K 



2^3 



-2Ka 



2k 



(87) 
(88) 



Z? = + 



1 /a^ 3a2 3 \ ,1 /2a3 3a \ 
-'x-^ + ^J«-2Aa+-^^-^Jcos2Aa 



1 /^&^_ 36^_ 

A3 2A5 ) 

1 /a2 



A 

G = + 



1 /263 
sin 2A6 - - I — 
4 



3^ 

A2 2"A4 



cos 2A5 



. I X ^ I sin 2Aa + cos 2Aa 
4 V A 2A3 y 4A2 

\ (\? \ \ b 

T — TTq- 2A6 - — r cos 2A5 

4 V A 4A3y A2 

F = -;^(sin2Aa-sin2A6) 
4A 



(89) 



(90) 



(91) 



J = + 



1 fa 



3aV 3 \ „. 1 /2a3 



A A3 +2A^;^°^'^'^-21a^"A^'""'^" 



3a 



36' 



2 V A 



^3+^Jcos2A6+^ 



1 /263 36 



A2 A4 



sin2A6 



(92) 
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K = + - { —7T COS 2Aa — — it sm 2Aa 



2 \ A 2A3 / 2A2 



-2U"^J -^2A6+^sm2A6 

L = ^ (cos 2Aa - cos 2A6) (94) 

M = - + + ^ j smh2,6 - 2 + 2?J -^^2,6 (95) 

^=2(7 + ^; sinh 2g6 - ^ cosh 2g6 (96) 



17 



5 Discussion and Conclusion 



In summary, a new method, namely, the Extended Wronskian Determinant Approach 
has been suggested to calculate the energy spectra and the wave functions of ID Dirac 
equation. An integral equation which can be solved by iterative procedure has been 
established to study the inhomogeneous Dirac equation with perturbation potential. 
Using the solution of the corresponding homogeneous Dirac equation as the zeroth wave 
function, we can calculate the wave function and energy spectrum of the perturbative 
Dirac equation up to any order in principle. Tow kinds of exact solutions including 
solution of one square-well potential and double-square-well potential have been em- 
ployed as the zeroth approximations to study the onc-wcll potential and double-well 
potential respectively. For one-well case, the first order iterative approximation to the 
energy levels and the wave functions has been obtained. For double-well case, we have 
got the energy level approximation up to the first order. These two examples confirm 
that the EWDA is a successful method to investigate the ID Dirac equation. 

In fact, it is not necessary to adopt the exact solutions of homogeneous ID 
Dirac equation as the zeroth order wave functions. Instead of exact solution, we 
ca use a trial wave function. And we can tisc the variational method for the non- 
perturbative Hamiltonian to determine the best choice of the parameters in the trial 
wave function first, and then use the trial wave function with the best parameters as 
the zeroth order wave function to study the perturbative Dirac equation by EWDA. 
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Figure Captions 

Figure 1: one- well potential 

Figure 2: perturbed one- well potential 

Figure 3: square-double-well and double-well potentials 
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